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' Abstract 

The muon flux at the South-Pole was measured for five zenith angles, 0°, 15°, 35°, 
82.13° and 85.15° with a scintillator muon telescope incorporating ice Cherenkov 
tank detectors as the absorber. We compare the measurements with other data and 
^ ! with calculations. 
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c/3 . 1 Introduction and motivation 

> 

^ ■ The IceCube neutrino telescope [1] under construction at the South Pole in- 

■ eludes a surface array, IceTop, for calibration and tagging of cosmic-ray in- 

duced background. Together, the surface array and neutrino telescope con- 
stitute a three-dimensional air shower array that will also be used to study 
the cosmic-ray spectrum from 300 TeV to 1 EeV [2]. IceTop consists of ice 
Cherenkov detectors with two tanks per station near the top of each IceCube 
string. The surface detector, when complete, will form a kilometer-square air 
shower array with nearest neighbor spacing between stations of approximately 
125 m. The tanks contain clear ice with a surface area of 2.4 m^ and depth 
of 0.9 m, and each is instrumented with two IceCube digital optical modules 
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(DOMs) partially embedded in the ice and facing down. Four two-tank sta- 
tions with a total of 16 DOMs have been taking data since they were deployed 
in the 2004-2005 austral summer season [3]. 

Water Cherenkov detectors were used as the primary detector elements in the 
Haverah Park surface array [4] and they make up the ground array of the 
Auger Project in Argentina [5]. Following an early test at the South Pole [6], 
development of a frozen Cherenkov detector as an element of an air shower ar- 
ray was undertaken in connection with the South Pole Air Shower Experiment 
(SPASE) [7]. Use of ice tanks as detectors for IceTop requires understanding 
the background radiation and in particular the muon fluxes at the South Pole. 
As in the Auger experiment [5], the pulse-charge distribution of single tank 
hits with its characteristic muon peak will be used for detector calibration 
and monitoring. In normal data taking muons will not be tagged by a muon 
telescope, so their fluxes must be well understood as a function of zenith angle 
on-site. In addition, muons are one of the major components of air showers, 
becoming especially important for showers at large zenith angle. Therefore 
detailed study of detector response to muons is needed to design the data 
acquisition system and interpret waveforms generated by air showers. 

Measurements of the muon flux are also of intrinsic interest, for example, in 
connection with understanding the production spectra of atmospheric neu- 
trinos. Although many muon flux measurements have been made throughout 
history (see the review by P.K.F Grieder [8]), most of them were done for 
vertical or nearly vertical muons. Only a few experiments measured nearly 
horizontal muons above a few GeV, for example at 3220 m a.s.l [9] and 3250 
m a.s.l [10,11]. In this work, using the test tanks deployed at the South-Pole 
as the absorber, we made a controlled and precise measurement of the flux 
of muons with minimum kinetic energy of several hundred MeV to about one 
GeV and zenith angle from vertical to nearly horizontal at the South Pole 
(~ 700 g/cm?, or ~ 3100 m a.s.l.). The details of the setup are given in §2 
and the analysis and results of our measurements in §3. To give a systematic 
view of the muon flux at this altitude, some data from other measurements 
are collected in §4. We conclude in §5 with a comparison to calculations. 



2 Details of the setup 

Most muons at the surface of the Earth are produced high in the atmosphere 
from the decay of charged pions generated in the interactions of high energy 
cosmic rays with atmosphere nuclei. The muon flux varies with altitude and 
has a strong zenith angle dependence, reflecting the convolution of their pro- 
duction spectrum with their energy loss and decay. At large zenith angles 
the muon flux becomes very small, as most of the muons generated in the 



2 



atmosphere decay before reaching the detector. As a consequence, the exper- 
imental setups for measuring the nearly vertical and nearly horizontal muon 
fluxes are different. Near the vertical a simple coincidence requirement is suffi- 
cient, whereas the measurement at large zenith angles requires more complex 
trigger and analysis procedures because of the low event rate. 

We show in Fig. 1 the configuration used for measuring the flux of nearly 
vertical muons. From December 10 to 24, 2000, the integral muon flux was 
measured at zenith angles of 0°, 15° and 30°. The absorber was a small ice 
Cherenkov detector installed at the South Pole inside the boundary of the 
SPASE-2 array [12]. It is a cylindrical polyethylene tank of area of 1.14 and 
height of 1.24 m. The inside of the tank is lined with white, diffusely-reflecting 
Tyvek (type 1025D). After flUing with the South-Pole station drinking water, 
two analog AMANDA optical modules (OMs) [13], with a separation of 0.51 m, 
were mounted facing down symmetrically off-center with their photo cathode 
region completely submerged. After the tank was frozen, the resulting ice 
depth was 0.99 m. A muon telescope consisting of three 0.2 scintillators, 
two stacked on top of the tank and one underneath, identifies penetrating 
muons by requiring a coincidence among SI, S2 and S3 within 50 ns. Electric 
pulses from the two OMs and the three scintillators were transmitted through 
~ 100 m of twisted pair and RG-8 coaxial cables, respectively, to the central 
SPASE-2 building that houses the electronics and data acquisition system. 
Three scintillator signals first enter a Phillips 711 discriminator, then go to 
a Phillips 755 logic unit. Both the singles rate and the coincidence rate were 
measured by a Jorway 1880B scaler. A 1 GHz Tektronix digital oscilloscope 
and a Linux PC were used to digitize and read out the waveforms from the 
two OMs through PCI-GPIB. 
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Fig. 1. A sketch of the ice Cherenkov detector (as the absorber), the muon telescope 
scintillators and the electronics in the flux measurement of 0°, 15° and 30° muons. 
Three fold coincidence was used to effectively eliminate the accidental coincidence 
rate and other background triggers that may form a two fold coincidence by chance. 



-J 



icriminator coincidence 
units 



(rq 



Linux PC 



scaler 



3 



The configuration used for the measurement of horizontal muons is shown in 
Fig. 2. Because of the low rate of horizontal muons the coincidence rate of 
two scintillators forming a simple muon telescope is dominated by triggers of 
local small air showers. Reducing the background requires special attention 
to the detector setup and to the data taking. First, as shown in the setup 
sketched in Fig. 2, two scintillators marked 'veto-Sl' and 'veto-S2' were used 
in anti-coincidence with SI and S2 respectively. Second, in addition to the 
waveforms from the two OMs in the tank, the trigger pulse and the traces 
of the discriminated pulses from the two muon telescope scintillators (SI and 
S2 on the figure) were recorded by two digital oscilloscopes to provide the 
relative timing information among these pulses. The cuts used to filter the 
data are described in the next section. The singles rate of each scintillator 
and the coincidence rate with and without anti-coincidence were monitored 
by a scaler that has a neghgible dead time at the working rates. The rates 
measured by the scalars were then used to correct for the dead time in the 
signal digitizing and reading system. 
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Fig. 2. A sketch of the detector layout (top) and the electronics (bottom) in the 
flux measurement of 82.13° and 85.15° muons. Scintillators 51 and 52 were used to 
form a muon telescope with veto-Sl and veto-S2 in anti-coincidence. 



The data for muon fiux measurement at zenith angle 82.13° and 85.15° were 
taken from December 13 to 21 in 2002, during which some time was spent 
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to rearrange the cables and detectors in the snow and to survey the detector 
positions. In this measurement, the tank was a larger cylindrical test tank 
(3.94 X 1.23 m) installed inside the boundary of the SPASE-2 array in 
2001. The preparation and installation of this tank were similar to the smaller 
tank. The two OMs are positioned symmetrically about the center with a 
separation of 1.07 m. The ice depth in the tank is 1.06 m. Fig. 3 shows a plan 
view of the coordinates of the tank and scintillators within the SPASE array. 
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Fig. 3. The coordinates of the detectors in the nearly horizontal muon measurement. 
S3, S6 and S7 are three of the SPASE-2 stations. Coordinate (X,Y) in meters are 
given in the SPASE-2 reference frame. 



3 Data analysis and muon flux at the South Pole 



3.1 Nearly vertical muons 



The flux of penetrating muons at three near vertical zenith angles is sum- 
marized in Table 1. Only the statistical errors are given. Column three gives 
the acceptances of the muon telescope at these three zenith angles. The ac- 
ceptances were obtained by a Monte Carlo simulation that requires 5.0 x 10^ 
particles to pass through all three scintillators for each case. The last two 
columns show respectively the integral fluxes and the minimum kinetic energy 
(and momentum) required for muons to pass through all the three scintillators. 
The flux is obtained as the measured rate of events divided by the acceptance. 
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Table 1 

Rates and fluxes of muons near the vertical at the South Pole. The muon thresh- 
old kinetic energy was calculated from the table of the continuous-slowing-down- 
approximation (CSDA) muon range in water [14]. The ice density used in the cal- 
culation is 0.92 g-cm~^. 



zenith(min,max) 


rate 


acceptance 


integrated flux 


Emin {Pmin ) 


n 


{Hz) 


(cm^ • sr) 


(cm^ • s ■ sr)~^ 


MeV{MeV/c) 


0.0(0, 16.5) 


2.22 ±0.01 


1.26 X 10^ 


(1.76 ±0.01) X 10-2 


246.0(335.4) 


15.0(0.4, 28.6) 


1.80 ± 0.01 


1.07 X 10^ 


(1.68 ±0.01) X 10-2 


263.5(353.7) 


35.0(25.8.44.2) 


0.72 ±0.01 


1.81 X IQi 


(1.49 ± 0.02) X 10"- 


311.7(403.8) 



3.2 Nearly horizontal muons 



Because of the low horizontal muon flux and the relatively high rate of local 
small air showers, the coincidence rate between scintillator 5*1 and S2 on 
Fig. 2 is dominated by background even when the two veto scintillators are 
included in the trigger. Figures 4 and 5 show the distributions of time of flight 
(TOF) as recorded by scintillators SI and S2, {TOF = Tgi - Ts2), with 
no cuts (the top histogram) and after various cuts have been applied to the 
data. The physical distance from 5*1 to 5*2 corresponds to a time differences 
of Tsi — Ts2 = —30.5 ns for the setup at 82.1° and to —63.5 ns for the larger 
angle. In the top histogram (without cuts) there is little or no sign of a peak at 
the expected times. Therefore, simply counting all events on a TOF window 
would greatly overestimate the true horizontal muon flux. The following series 
of cuts was applied to reduce the background: 

• cuti. For each triggered event, scanning the waveforms recorded by the 
digital oscilloscopes and choosing those with exactly one hit in the two 
triggering scintillators and two OMs in a time window of 500 ns. 

• cut2'- Choosing events in which the peaks of the pulses in the two OMs in 
the tank are close in time abs{ToMi —T0M2) < 14 ns. The 14 ns is related to 
the slow rise time and the signal fluctuations in the tank. This cut selects 
those events of which the Cherenkov photons generated in the tank give 
rise to the signals in the two OMs evenly and simultaneously. The first two 
cuts arc both designed to reduce the number of events caused by showers 
with several particles spread in time. In addition to single particle events, 
however, events with a few particles arriving within a time window less than 
the detector resolving times can still pass these cuts. 

• cwta: Choosing events in which abs{ToMi — Ts2 — 234 ns) < 16 ns. Here, T52 
is the time of the discriminated pulse from the trigger scintillator 5*2. The 
average cable delay between the OMl signal and the signal from scintillator 
S2 recorded on the oscilloscope is 234 ns. It was determined by a system 
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calibration with vertical muons in which we used the same cable, the same 
trigger threshold, etc. The 16 ns is one half of the rise time of the OM pulse, 
the fluctuation in the tank signal as seen in OMl and the fluctuation in the 
signal from scintillator S2. This cut gets rid of those triggers in which the 
pulse generated in 5*2, which is located next to the tank, is too far in time 
from the signal in the tank. 
• cut4: Choosing events in which abs{Tsi + Ttof + Tdeiay) < ^if)- Here, 
Tsi is the time (on the oscilloscope trace) of the discriminated pulse from 
scintillator SI. In the setup, it is a negative number (i.e. to the left of the 
trigger moment on the oscilloscope trace) for single muon event. Ttof is the 
muon time-of-fiight from SI to S2, which is 30.5 ns and 63.5 ns for muons of 
6 = 82.13° and 9 = 85.15° in this setup. Tdeiay is the 50 ns trigger delay set 
in the digital oscilloscope. For single muon events, abs(Tsi + Ttof + Tdeiay) 
is the time jitter relative to the DAQ trigger moment which is ns on the 
oscilloscope trace. In the analysis, 6{t) was chosen 4.5 ns, which corresponds 
to the time jitter in the scintillator and the electronics. 
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Fig. 4. The TOF {Tsi — Ts2) spectra of triggered events in ^ = 82.13° muon mea- 
surement. The cuts used for selecting the events under each histogram are marked. 
See more details about the cuts in the text. 



7 



85.15° 



1 -v 



10' 




-25 25 



75 100 
Tsi - (ns) 



10 - 



cut4 only 



n 



all four cuts 




-100 -75 -50 -25 25 



50 75 100 

Tsi - (ns) 



Fig. 5. The TOF spectrum of triggered events in 6 = 85.15° muons measurement. 

In Fig. 4 and 5, from top to bottom, the 51 to 5*2 TOF spectrum are shown 
for triggered events without any cut, and events left after the apphcation of 
cuts 1 to 4. The following features on the TOF spectra are noted: 

(a) The events left after the application of the first three cuts start to show a 
clear peak at the correct TOF position on the spectrum although none of 
these cuts makes use of the time-of- flight from 5*1 to 5*2; 

(b) Although most of the events that survive the cut^ (and no other cuts) dis- 
tribute themselves around the expected muon flight time on the TOF spec- 
trum, several of them still fall outside the expected time window in Fig. 4, 
and more in Fig. 5; 

(c) All the events left after all four cuts show themselves at the expected posi- 
tion as on an ideal TOF spectrum of single muon events; 

(d) The events in the peak with all four cuts applied correspond to those events 
that make up the peak on the distribution with cuts 1, 2 and 3. 



We therefore conclude that the events left after application of all four cuts are 
true horizontal muon events. It is interesting to compare their waveforms with 
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vertical muon waveforms as seen by the two OMs in the same tank. We note 
two features in the data: 

(1) The average amphtude of these pulses is about 2 times as big as that of 
the vertical muons. 

(2) The average rise time of the these pulses in the tank is 3ns longer than 
that of vertical muon pulses. The fall time of the average waveform is 
nearly the same as that of the average waveform of vertical muons. 

A GEANT-4 simulation of the tank response was carried out for both vertical 
muons and muons of zenith angle 82.13° and 85.15° [15]. The two features 
listed above are well reproduced in the simulation. They are directly related 
to the longer track length of the near horizontal muons in the tank. 

Taking the events left after all four cuts as the single muon events, we sum- 
marize the measurement in Table 2. 



Table 2 

Summary of the horizontal muon measurement at the South-Pole. The data taking 
time excludes dead time. The systematic errors are from the uncertainty in the 
cut parameters. They were obtained by a systematic comparison of the SI to S2 
time difference spectra under four cuts with different parameter values. The muon 
threshold kinetic energy was calculated from the table of the continuous- slowing- 
down-approximation (CSDA) muon range in water [14]. The ice and the snow 
densities used in the calculation are 0.92 g-cm~^ and 0.4 g-cm~^. The latter has an 
error of ±5%. 



zenith coverage (°) 
time-of-flight (ns) 
number of events 
acceptance (cm^ • sr) 
data taking time(sec) 
EminiMeV) {p„nn{MeV/c)) 
statistical errors 
systematical errors 
flux (cm^ • s ■ sr)"^ 



82.13 ±2.97 
30.5 
173 
4.6881 
2.1227 X 10^ 
899 ± 45(999 ± 50) 
±7.6% 

+14% 
-13% 



85.15 ± 1.44 
63.5 
33 
1.0944 
2.5600 X 10^ 
1752 ± 87(1855 ± 93) 
±17.4% 

+3% 
-12% 



(1.74l°;24 ± 0.13) X 10-4 (l-18lo;?4 ± 0.21) x 10"^ 
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4 Muon fluxes by some other experiments at an altitude close to 
the South-Pole 



There are several measurements of the vertical muon flux at high altitude 
that are performed at different muon threshold energies, geomagnetic lo- 
cations and solar epochs. At 3220 m a.s.l, the vertical intensity of muons 
above 2 GeV was measured with GM hodoscope by Shen and Chiang [9]. 
The reported flux is (4.9 ± 0.2) x 10~^(cm^ ■ s ■ sr)~^. For muons of energy 
> 162.9MeV, the vertical integral intensity at the South-Pole altitude is 
(2.61 ± 0.01) X 10-2 (cm^ ■ s ■ sr)'^ according to the curve on figure 2.130 
in reference [8], which was made from the measurement by Blokh et al. [16]. 
Vertical muons at 2960m a.s.l. and 3250m a.s.l. were measured by AUkofer 
and Triimper [17] and Kocharian et al. [18]. By integrating the differential 
momentum spectrum on figure 2.137 in reference [8], the integral muon inten- 
sities were found to be /2960m(> 203 MeV/c) = 1.16 x 10~'^(cm'^ ■ s ■ sr)~^ and 
l3250m{> 4:^4: MeV/c) = 1.58 x l(}~''^{cm'^ ■ s ■ sr)^^. Some of the differences 
between these results can be explained by the different altitudes and muon 
threshold energies. The difference in altitude between Ref. [17] and South Pole 
corresponds to about 30 g/cm^, which increases the energy threshold by 60 
MeV. In addition, there is also a decrease of the muon flux because of muon 
decay. About 15% of the muons at threshold decay in 300 m. These effects, as 
well as the solar epochs, do not explain all the differences among the quoted 
measurements. Our result is in the middle of the scattered points obtained by 
those previous work and agrees the best with the result of Ref. [18]. 

Few experiments have been done for low energy nearly horizontal muons with 
a small zenith bin. For an angular interval from 78.4° to 90°, corresponding 
to a mean of 86.2°, the intensity at 3220 m a.s.l was reported to be [9]: //i(> 
2 GeV) = (9.68 ± 0.32) x 10-^{cm^ ■ s ■ sr)'^. Because of the large angular 
range and the steep dependence of the muon flux on the cosine of the zenith 
angle this number which is lower than the current measurement is difficult to 
interpret. 

Asatiani et al. [10] carried out a spectral measurement at 3250 m a.s.l with 
a large magnetic spectrometer. The muon differential spectrum covers the 
range from 10 GeV/c to 2000 GeV/c. The muon incident angle is from 80° to 
90° with the mean zenith angle 84°. The presented spectrum was normalized 
to the calculation by Ashton et al. [19] at 50 GeV for the zenith angle 84°. 
Integrating over the differential energy spectrum in figure 2.139 in reference [8] 
gives an approximate intensity of (8.85 ± 0.02) x 10~^ {cm'^ ■ s ■ sr)~^. 

With the same instrument, Asatiani et al. [11] reported the differential muon 
spectra at narrow zenith angle intervals covering from 80° to 88°. The integral 
muon fiux in units of (cm^ • s ■ sr)~^ above a given energy threshold for four 



10 



zenith bins are: 80° - 82°, />i6.5Gey = 1-27 x 10"^; 82° - 84°, I>ie.5GeV = 
1.01x10"^; 84°-86°, I>25.8GeV = 5.60x10-^ 86°-88°, I>u.8GeV = 4.40x10"^ 

Because of the finer binning in the zenith angle we can correct these numbers 
for the difference in muon threshold and compare to the current measurement. 
The correction was done using the calculation of Ref. [20] for the fluxes of 
muons at cos^^ = 0.25, 0.15 and 0.05 at the sea level. The muon fluxes from 
this calculation were interpolated for the average cos 6 and energy threshold in 
the measurement and compared to the fiuxes above 2 GeV energy threshold. 
The scaling factors for the four points were found to be 1.84, 1.60, 1.77 and 
1.22 respectively. The data of Ref. [11], after being scaled to muon energy 
threshold of 2 GeV, are shown with open squares in Fig. 6. The agreement is 
excellent. 



5 Compcirison with simulation 

The measurements were compared with simulations. For the vertical and near 
vertical muons, a Monte Carlo simulation was carried out to calculate the 
fiuxes expected for the average atmosphere depth of 702 g/cm? during the 
muon measurement. The simulation used TARGET2.1 [21] and the cosmic 
ray flux of reference [20] for the epoch of solar maximum. The muon fluxes 
of 1.80 X 10~^, 1.73 X 10"^ and 1.68 x 10"^ {cm? ■ s ■ srY^ were predicted 
for the three angles. The measured muon fiuxes at 0° and 15° are in excellent 
agreement with the prediction, while at 35° the fiux is lower by 12%. 

The fiuxes of nearly horizontal muons by this work are compared to the the 
simulation of Ref. [20] . The current measurements are shown with solid squares 
in Fig. 6. The predicted angular dependence of muons with energy above 1 
(1.86) GeV from Ref. [20] are shown with a dashed (solid) line. One can see 
that the expected fluxes for the two muon energy thresholds at cos 9 less than 
0.2 become the same. This is because the horizontal muon energy spectrum 
has very few muons in this energy range. The good agreement between the 
measurements and the calculations also confirms that accounting for difference 
in the altitude of the measurements is not important for nearly horizontal 
muons. 



6 Summciry 

Using the prototype IceTop ice Cherenkov detector as the absorber, data has 
been collected for muons in five directions from vertical and nearly vertical to 
nearly horizontal. To select the rare nearly horizontal muon events, four cuts 
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Fig. 6. The measured integral flux of muons above the energy thresholds specified 
in Tables 1 and 2 are shown with full squares versus the cosine of the zenith angle. 
The shaded boxes around data points show the errors in the fluxes and the zenith 
angular coverage in the measurement. The open circles show the calculation for 
vertical and near vertical muons at the South Pole. The lines show the angular 
dependence of the predicted integral fluxes above the two energy thresholds. The 
open squares are the measurements of Ref. [11] after being scaled to a threshold of 
2 GeV. 

were developed in the data analysis. The procedure for selecting horizontal 
muons was crosschecked by a GEANT-4 simulation of the muon waveforms 
seen in the OMs. 

The integral muon flux of muons at 0° and 15° agrees very well with the simu- 
lation, which accounts for the atmospheric depth and the specific atmospheric 
density profile above the South Pole. The simulation of Ref. [20] is not done 
for the South Pole location, however this is not important for the high energy 
muons that survive at very large zenith angles. 
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